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Abstract 

We consider the phase transition in the system of n simuhaneously developing random walks on 
the halfline a; > 0. All walks are independent on each others in all points except the origin a; = 0, 
where the point well is located. The well depth depends on the number of particles (i.e. steps of 
different random walks) simultaneously staying at the point x — 0. We consider the limit n 3> 1 
and show that if the well depth grows faster than |nlnrt with n, then all random walks become 
localized simultaneously at the origin. In the conclusion we discuss the connection of the above 
defined problem with the phase transition in the copolymer chain with quenched random sequence 
of monomers considered in the frameworks of replica approach. 
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Introduction 



Problems dealing with the localization of random walks in potentials of various geometries compel 
much attention of different groups of scientists: both chemists and physicists. For the first group 
these problems are naturally connected with creating a new materials possessing specific technological 
properties (for instance, catalysts Q or ion containing and surface active substances [^). Meanwhile, 
for physicists the consideration of random walks localization remains since 70's the firing ground for 
elaborating the new methods of the investigation of a polymer absorption in different geometries, 
wetting phenomena and kinetics of chemical reactions. 

The phenomenon of polymer chain absorption is rather well understood at present both for homo- 
[^] and heteropolymer [^-[0] cases. The simple diffusive approach 0, |l^] provides com- 
plete understanding of an absorption of homopolymer chains as well as of block-copolymers in com- 
plicated geometries. These results can be beautifully explained also by scaling estimates |l^ . More 
advanced renormalization group methods |^] and power series analysis |jl^ applied to random chains 
with a disordered sequence of links are also widely used and give an exhaustive information about the 
thermodynamic properties of ideal polymers near the point of phase transition from the delocalized 
(Gaussian) to the localized (absorbed) fluctuation regimes. The similar mathematical formalism has 
been applied to a description of a phase transition in solid-on-solid (SOS) models with quenched 
impurities jlT]]. Many conclusions obtained by RG analysis in [§] correlate with the results in p 



Problems dealing with 2D wetting subject to a layer of periodic potential could be attributed to 
this theme as well. In [18| using a generalized transfer matrix method, we were able to find the exact 



solution for the critical depinning phase transition. 
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Finally, it is worthwhile to mention that the model considered below has some features of so- 
called "zero's range process" emerged in connection with traffic models. Also, similar set 
of problems appeared in course of investigation of the " autolocalization" phenomena |21] and in the 
theory of multicomponent chemical reactions on traps immersed in an immovable "matrix" ||2^, p3| . 

The present paper is organised as follows. In Section |^ we formulate the simple model of "catalytic 
absorption" and put forward the question of our main interest; Section Q is devoted to the solution 
of one- and two-particle problems; the consideration of the phase transition in the n-particle model 
for n > 2 is given in Section ^; while in Section Q we present conjectures about the connection of the 
multiparticle catalytic absorption with the random copolymer absorption at the point well. 



1 Multiparticle model of catalytic absorption at the point well 

Consider the one-dimensional lattice 2"^, i.e. the set of integer points on the halfline x > 0, {x} = 
{0, 1, 2, 3, . . .}. Take n independent A^-step simplest random walks on^"*", simultaneously starting at 
the point x = 0. All these walks are free, i.e. they do not interact neither with any external potential, 
nor with each others in all points except the point x = 0. At the point x = the interaction of a given 
random walk with the well is described by the potential U(xi,X2, ■ ■ ■ ,Xn) (the index i E lables 
the different random walks): 



U {xi,X2, ...,Xn) 



if Xi / for all i G [l,n] 

?7(i) if Xi = 0; Xj ^ for ah {j ^ i; j e [1, n]} 

if x^ = 0, Xj =0; Xfc / for all {k ^ j ^ i; ke [1, n]} (1.1) 

_ C/(") if xi = 0, X2 = 0, . . . , = 

< C/(^) < C/(2) < • • • < (1.2) 

^«^/(i); iG[l,n] (1.3) 



where 
and 



The equations ( |l.l| )-( pr3D reflect the fact that the depth of the potential well depends on the 
number of particles (i.e. steps of different random walks) simultaneously located at the point x = — 
see Fig|l|. Let us stress that U^^^ is the "individual" potential, i.e. the potential per one particle in 
a cluster of i particles simultaneously located at the point x = 0. The "integral" well depth for the 
cluster of i particles we denote as 

U{i) =iU^^ = if{i)- 

It is very straightforward to derive the recursion relations describing the process under investiga- 
tion. Consider the partition function (the number of trajectories), ZN{xi), of some arbitrary A^-step 
random walk on the halfline Xi £ [0, oo) starting at the point x = 1 and ending at the point Xi 
{xi gZ+). The function Z]\f{xi) satisfies the equations: 

ZN+liXi) = Z^iXi + ZN{Xi - 1) {Xi > 1) 

Z;v+i(x.) = e^(^i'-'^--'^")Z^(xi + l) (x, = 0) (1.4) 

ZN=o{Xi) = Sx^fl 
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The same "equations of motion" (as Eq. (|1.4| )) should be written for all functions Zi^{xi) where 
i =G [1, n] and the potential U (xi, . . . , Xj, . . . , x„) is defined in Eq.( |l.l| ). 

We are interested in the situation when the number of particles jumping onZ"*" (i.e. the number 
of different random walks) is very large (n 1) and the number of steps of each random walk tends 
to infinity [N — > oo). 

We expect that in the thermodynamic limit — > oo and for n S> 1 the absorption transition 
into localized state at the point x = is sensitive to the shape of the function f{n). Namely, if f{n) 
grows rather slowly with n (for the precise criteria of the growth see the Section ^) we could expect 
that the phase transitions occur in each individual random walk independently; while if f{n) is rather 
sharp function of n, various types of collective localization are expected. The interplay between the 
entropy loss and the energy gain in the localized states makes the phase behavior of the system under 
consideration extremely rich. 

The most attention in our work is paid to exact determination of the critical shape of the function 
/(n) at which the localization transition occurs simultaneously in all n (n ^ oo) random walks. 

To avoid further misunderstandings let us define precisely what we mean under the "localization" 
of the random walk. Let F{N) be the free energy of the iV-step random walk on the halflineZ^ with 
the point well of depth U at the origin x = 0. This model can be described via the diffusion-type 
equation with specific boundary conditions (compare to Eq. (|1.4| )) and there exists some critical value 
U = UcT > which separates two different behaviors of the free energy: 

1 , , f for < C/cr 

lim —F(N) = I 

^ [ A > for ?7 > C/cr 

The value Ucr we call the "localization transition point". It signals the separation of the highest 
eigenvalue A from the continous spectrum of the diffusion-type equation. 



2 Exact solutions of the catalytic absorption model for n = {1, 2} 



2.1 Solution for n = 1 

Let us start with the simplest situation when n = 1, i.e. we have a single simplest random walk on 
Z"^ interacting with the potential well at the origin. 



In this case Eqs.(l.l)^ (|l.4D degenerate to the following ones: 

U{x) = 



if X / 
?7(i) if X = 



(2.1) 



and 



Zn+i{x) = Zn{x + I) + Zn{x - 1) 
Z^+i(x) =e^(^)Zjv(x + l) 
Zn{x) = 
Zn=q{x) = 5xfi 



(x > 1) 
(x = 0) 
(x < 0) 



(2.2) 



Thus, for n = 1 we have only one particle jumping on the halfline and the problem is nothing else as 
a simple absorption at the "point well" located at the origin x = 0. 



The solution of the one-particle problem is very straightforward. Let us make the shift x ^ x + 1. 
Now the equations ( |2.2| ) can be rewritten as: 

ZN+iix) = ZNix-l) + ZNix + l)+5,,i(^e^^'^ -l)ZNix + l) (x > 1) (2.3) 

completed by the boundary Z]\f{x = 0) = and initial Zj\f=Q{x) = 6x,i conditions. 
Perform the substitution 

Zn{x) = 2'' Zn{x) 

and use the Fourier-Laplace transform 



(2.4) 



oo oo 



N=0 x=0 

irq 



Introducing the new variable k = —, we arrive at the following integral equation for the function 
ZNik,s): 



1 



1 



— Z{k,s) sin k = cos k Z{k, s) + sin k (e^^ ^ — 1 ) — / dk sin 2kZ(k, s) 

s s V ^ vr Jo 



1 



(2.5) 



After simple algebra we get 



Z{x = 2,s) = — / dksm2k Z{k,s) 
vr Jo 



vr J 


/ dk 
'o 


sin k sin 2k 
1 — s cos k 




1 - . (e^<^^ 


-0 


1 r „ sin 
vr Jo 1 - 


k sm2k 
- s cos k 



(2.6) 



The divergence of the function — / dk sm2k Z{k, s) occurs at the critical value U^^ = Ucv (corre- 

vr Jo 

sponding to the localization transition point) when the denominator in Eg. ( p. 61) is set to zero: 



(1) 



1 



1 sink sm2k 

^~ / T 

vr Jo 1 — s cos k 



(2.7) 



Thus, in the thermodynamic limit — > oo (i.e. when s — > 1) we have C/ir^ = In 2. 



2.2 Solution for n = 2 
2.2.1 General ansatz 

Consider now two particles jumping on the halfline, i.e. two trajectories of random walks simulta- 
neously developing in "time" t {0 < t < N). These trajectories are independent on each others in 
all points x = {1, 2, . . .} except the point x = 0. Let us call the current coordinates of 1st and 2nd 
random walks xi{t) and X2{t) correspondingly. The potential (Eq. (|l.lD ) reads now: 

{0 if xi / and X2 / 
[/(I) if {xi = and xs 7^ 0} or {xi / and X2 = 0} (2.8) 
if xi = and X2 = 
4 



Thus we have the 2D random walk in the first quarter of the (xi, a;2)-plane with specific boundary 
conditions. In order to simplify the corresponding equations, it is very convenient to choose the 
embedding lattice as it is shown in FigJ^ The equations for the joint partition function, Z]\f{xi,X2), 
in 2D-case (for n = 2) become more tricky than in ID-case (for n = 1). 

Performing the shift xi ^ xi + 1 and X2 — > a:2 + 1 (see Fig.^), we can derive the recursion relation 
for the partition function Z]\f{xi,X2) which is the two-dimensional extension of Eq.( |2.3| ): 

Zn+i{xi,X2) = Axj^,x2Zn{xi,X2) + 

6x1,1 [e^'^^^ - 1^ A^i^Nixi + 1,X2) + 

5x2,1 - Ij AxiZAr(xi,X2 + 1) + 

5xi,i<5x2,i (e^^*'^ - 2e^*'' + l) Zj,(xi + 1,X2 + 1) 

where 

Ax^^'(xi) = ^'(xi - 1) + ^'(xi + 1) (« = {!, 2}) (2.10) 

and 

Axi,X2*(a;i, X2) = ^(Xi - 1, X2 - 1) + ^(Xi - 1, X2 + 1) + *(xi + 1, X2 - 1) + ^(Xi + 1, X2 + 1) (2.11) 



The equation ( |2.9D is valid for xi > 1, X2 > 1 and should be completed by the boundary and 
initial conditions: 

Zfqixx = 0, X2 > 0) = Zn(xx > 0, X2 = 0) = Zn{xx = 0, X2 = 0) = ^ 
Zjq=Q{xx,X2) = 4i, 1(5x2,1 

It is easy to check that Eqs.( |2.9| )-( |2.12| ) describe properly our model (taking into account the shift 
3^1,2 Xi,2 + 1): 

(i) For xi > 2 and X2 > 2 we have a simple diffusion without any interactions; 

(ii) For xi = 1 and X2 > 2 (or for xi > 2 and X2 = 1) the Boltzman weight P^^^ of the random 
walk steps located at the point xi = 1 (or X2 = 1) reads 

/?« = (e^'''-l)+l = e^*''; (2.13) 

(iii) For xi = 1 and X2 = 1 the Boltzman weight P^"^^ of the steps of two different random walk 
simultaneously located at the point (xi = 1,X2 = 1) reads 

^(2) ^ (^2C/(^) _ 2,^w + i) + 2 (e^'^' - l) + 1 = (e^'^')' . (2.14) 

Let us search the solution of Eqs.(p.S|)-( |2.12| ) using the following ansatz 

^iv(xi,X2) = Zjv(xi)Z^^'=(x2) + Z^(x2)Z^^''(xi) + PF^(xi,X2) (2.15) 

where: 

a) The "one-particle" partition function Z]^{xi) (i = {1,2}) is the same as in Eq.( p.3D ; 
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b) The term Z^^'^'^{x) is defined as follows 

Z^-(x) = A,. Zfr-(x) (x > 1) 

Z^'^^x) = (x = 0) (2.16) 

and has the sense of "free" partition function without any potential on the halfline x > 0; 

c) The last term W]\f{xi, X2) describes the "absorption at the main corner^' and has the "nonmulti- 
plicative nature" being determined by the recursion relations: 

WN+lixi,X2) = A^^^^^WNixi,X2) + Sx^,iJ2Wn{xi + 1,X2 + 1) (2.17) 

where 

' WNixi = l,X2>2)=WNixi>2,X2 = l)=WN{xi<0,X2<0)=0 

[Z.Lii) 

Wn=o{xi,X2) = 6x1,1 5x2,1 

and the "main corner" Boltzmann weight 

^2 = e'U''' - 2e^''' + 1 (2.19) 

is chosen such that Eq.( 2.15 ) reproduces the right weight (3^"^^ at the point {xi = 0,X2 = 0) (see 
Eqs.dJ) and (|1|)). 

The equations ( |2.17] )-( [2.18| ) correspond to the situation shown in Fig.^a, where we marked the 
point at which the Boltzmann weight 72 is located by the sign [•] . 

As it has been stated in Section |l| we are interested in the determination of the phase transition 
point of the function Z]\f{xi,X2) (for n = 2 and 1). Generalising the arguments of Section 2A it 
is easy to verify that the transition to the localized state occurs at the divergence point of the function 
Z{xi = 2,X2 = 2,s) (as in Eq.(U)): 

Zn^oo{xi = 2,X2 = 2) = \ I I dkidk2sm2ki sm2k9Z{ki,k2,s ^ 1) (2.20) 

vr^ Jo Jo 

where 

^7V(X1,X2) =4^Z7v(xi,X2); Wn{xi,X2)=4''Wn{xi,X2) (2.21) 

and 

00 00 00 

Z{ki,k2, s) = ^ ^ ^ sinxifci sinx2A;2 2^Ar(xi, X2) 

Af=0 x-i=0 3:2=0 

(compare to Eq.( |2.4| )). 

Substituting the ansatz ( p.l5| ) into ( p. 201 ) we see that the diverg ence of the function Z ^^{xi — 
2, X2 = 2) is determined: 

• either by the divergence of the "one-particle" function Zj\i^^{xi = 2) (i = {1, 2}) 

• or by the divergence of the "main corner" part W]\f-,ooixi = 2,X2 = 2) 

and everything depends on the fact which of these two functions diverges first (i.e. for smaller values 
of {C/(^) and C/^^)}) in the phase space of parameters {U^^\ U^"^^} (recall that according to Eq.( |l.2[) 
we have a restriction < U^^^ < C/*^^^). 



The behavior of the function ZN{xi) for i = {1,2} has been well studied in Section 2.1; so, let us 
concentrate the efforts on the solution of Eqs.( ^.17|) -( p.l8D . 



2.2.2 Phase transition point of the "main corner" part Wjv(xi,X2) 

The partition function W{N) = Wi\f{xi = 1,X2 = 1) where is even, it is possible to represent in the 
following form. Let us introduce the auxiliary functions: 



1. f2(iV) = flpf^xi = 1,X2 = 1) — the number of closed paths, starting at the point (xi = 1,X2 = 1), 
finishing at the same point after steps for the first time and satisfying the boundary conditions 
in the Fig.|a; 

2. V{N) = V]s[{xi = 2,X2 = 2) — the number of A^-step closed paths, starting and finishing at the 
point {xi = 2,X2 = 2) and satisfying the boundary conditions in the FigU^; 

The following identity can be derived immediately: 

n{N + 2) = V{N) = 4^V{N) (2.22) 



The function W{N) admits the representation 

N/2 k 

W{N) = J2{p(^^f J2 (2.23) 

k=l {Ni+...+Nk=N} 1=1 

where the upper limit in the first sum can be set to infinity because the condition A'^i + . . . + N/. = N 
ensures the right cut of the sum; /J^^^ = 72 + 1 is the Boltzmann weight of the potential well at the 
origin. 

Introducing the Kronecker (5-function: 

I r dz J x^O 
d[x) = : 



2TTi J z^+^ 1 1 x = 
where x = — A'^i — . . . — A^^^, we may rewrite Eq.( |2.23| ) as follows: 



> ' I, '\ f 00 00 fc 

W{N) = ^ id^z-^-i ^ ... ^ J]n(iV.) 

fc=0 Afi=2 Nk=2 i=l 

^^'"^ 1-/3(2) Y^z'^m) 

N=2 



The appearance of the pole in the last expression signals the separation of the localized mode from 
the continuous part of the spectrum corresponding to the function W{N). Using the identity ( ^.22 ) 
we get the equation on the transition point (for A^ — > cxd): 

= hm^ 00 ._^N = !™ (2-25) 



E l m) E-'^V(Ar) 

N=0 ^ ^ N=0 



where s = —. 

4 
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The function Vj\[{xi,X2) satisfies the recursion relation in absence of any potentials: 

VN+lixi,X2) = ^xi,X2 Vn{,Xi,X2) (2.26) 
with the boundary conditions shown in Fig.^b: 

Vn{xi = 1,X2>2) = Vn{xi >2,X2 = l) = Vn{xi = 1, X2 = 1) = 
Vn=o{xi,X2) = (5x1,2 (5x2,2 



Using the 2D-Fourier and Laplace transforms and performing the back shift in Eq.(2.26) {xi — > 
xi — 1, X2 — > X2 — 1} we arrive at the standard equation for the function VAr(A;i, k2, s): 

, N sinfci sinA;2 

VN{ki,k2, s) = ^ V 2.28 

1 — s cos ki cos k2 

Thus, we get the final expression: 



TT TV 



E s-ViN) = [dk, dk2 \ (2.29) 

TT^ J J I — S cos kl COS k2 



N=0 







Evaluating the last integral at the point s = 1 we obtain the following numerical value for the "main 
corner" transition point: 

= ^2 + 1 = - 3.65979 (2.30) 

4 — vr 

2.2.3 Phase diagram for 2 particle system 



Collecting the results of the Sections 2.2.1 and 2.2.2 we conclude that the localization transition from 



the delocalized to the absorbed state is determined: 



either by the equation /3cr = e = 2 (the "one-particle" part); 
part). 



or by the equation =72 + 1, i.e. e^^*^' — 2e^*^' + 2 = (the two-particle "main corner" 



The complete phase diagram is drawn in Fig.|5[ For U^^^ > ln2 the influence of the potential C/^^^ 
disappears, while for < f/^^-* < ln2 ~ 0.693 the "induced" localization can occur for 

t/(^)>ilnf2e^'%-^-2 
2 V 4-7r 

just due to the simultaneous interactions between two random walks at the origin. 

It is worthwhile to mention that the "collective" localization can appear even for U^^^ < 0. Suppose, 
for example, such situation: 

[/{!) = -00; > 

i.e. each particular random walk cannot visit the origin, but there is an energy gain for two random 
walks to stay at the origin simultaneously. In this case for 



[/{2)>[/(2) = iln('^_2 

2 V4-7r 



0.253 



the trajectories become localized. 



3 Phase transitions in catalytic absorption model for n > 2 



3.1 Solution for n = 3 



The equation for the 3-particle partition function Zn{xi,X2,X3) reads: 

Zn{xi,X2,X3) = Aj:^^x^^x,^Zn{xi,X2,X3) + 

Sxi,l (e^^'' - l) Ax2^.j;^Zn{xi + 1,X2,X3) + 
Sx2,l (e^*^^ - l) A^i,x3Zn{xi,X2 + 1,X3) + 
5x3,1 (e^*'' - l) Axi,X2Zn{xi,X2,X3 + 1) + 



<^Xl,l<Jx2,l 

5x2,1^x3,1 (e 

5x1, 1^X2, 1^X3,1 



,2C/(2) 
2C/(2) 



2e 
2e 



2e 



(7(1) 
[/(I) 



+ ij A^^Zn{xi + 1,X2 + 1,X3) + 
+ A^^Z^ixi + l,X2,X3 + 1) + 
+ l) A^^Zn{xi,X2 + 1,X3 + 1) + 



,3C/(3) 



(3.1) 



1 ZAr(xi + 1, X2 + 1, X3 + 1) 



Generahsing the ansatz (2.15) to the 3D-case we may search the solution of Eq.(|3.1|) in the form: 

ZNixi,X2,X3) = 2'Ar(xi)Z^^°(x2,X3) + Zn{x2)ZP"'{xi,X3) + Zat (x3)Z^^^(3;i , rE2) + 

Wn{xi,X2)Z%--{x3) + Wn{xi,X3)Z%^%X2) + Wn{x2,X3)Z%-'^(xi)+ (3.2) 
Wn{xi,X2,X3) 



where: 



a) The partition function W]\f{xi, Xj) {i ^ j; {i,j} € [1,3]) is the same as in Eq.( |2.17] ) 



b) The functions Zff^°{xi) and Z^^'^{xi,Xj) {i ^ j; {i, j} G [1,3]) describe the contributions of one- 
and two-dimensional free random walks with zero's boundary conditions (compare to Eq.( 2.15D ); 

c) The function iyAr(xi, X2, 2:3) is determined as follows: 



Wn +l{xi,X2,X3) 



iXi,X2,X3 



Wn{xi,X2,X3) + 5x-i,l 5x2,l5x3,l73^Af(2;i + 1,X2 + 1,X3 + 1) 



Wn{xi = l,X2> 2,X3 > 2) = WNixi >2,X2 = 1,X3 > 2) = WNixi > 2,X2 > 2,X3 = 1) = 
Wn{xi = l,X2 = 1,X3 > 2) = Wn{xi = l,X2> 2,X3 = 1) = Wn{xi > 2,3;2 = 1,X3 = 1) = 
^ WAr=o(xi,X2) = 5x1,1 5x2,1 

(3.3) 



with 73 = e 



3C/(3) 



1. 



In order to check the validity of ansatz (^^) let us compute the values of the Boltzmann weights 
at the boundaries: 

1. At the planes: 



Xl=0 {X2 >0,X3> 0) 

X2 = (xi>0,X3>0) \ : = (e^^' 



X3 = (xi > 0,X2 > 0) 



1+1 



(3.4) 
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2. At the lines: 



{xi=0,X2 = 0) X3>0 _ (2) (1) 



(xi=0,X3 = 0) X2 > 
(X2 = 0,X3 = 0) Xi > 



/3(2) = e2C/(^' _ 2e^'^' + 1+ 
(e^<^' - l) + 1 = e 



2(7(2) 



(3.5) 



3. At the "main corner" point: 



(X, = 0, X2 = 0, X3 = 0) : ^ _ ^ ^J,, (3.6) 



The arguments of the Section |2.2.2| can be easily extended to the 3D-case. Skipping the inter- 
mediate computations we bring the final result for the value of the Boltzmann weight at the "main 
corner" transition point: 

/j(3) =l±.J J Jdk,dk2 dh r^T^^lTl^^l 1 « 7.856 (3.7) 





The localization phase transition is now determined: 

• either by the equation (3ct^ = 2, i.e. 

U^P=ln2 (3.8) 

(the one-particle part); 

• or by the equation = 72 + 1 = 1 , i-e. 

4 — TT 

t/(2) = iln(2e^<%^-2) (3.9) 

(the two-particle "main corner" part); 

• or by the equation = 73 + 1 = 7.856, i.e. 

= i In (7.856 + 3e2^*'' - 3e^''') (3.10) 
(the three-particle "main corner" part). 

3.2 Solution for arbitrary n 

The general ansatz for arbitrary values of n reads 

Zn{xi,X2-, . . . ,Xn) = .^Ar(xi)Z^'''^(x2, . . . , x„) + all permutations+ 
TyAr(xi, X2).^^'^''(x3, . . . , x„) + all permutations+ 
WAr(xi, X2, X3)Z^°'^(x4, . . . , x„) + all permutations+ 
••• + 

Wn{xi, . . . , x„_i)Z^'^'^(x„) + all permutations+ 
Wn{xi, . . . ,x„) 

10 



(3.11) 



Eq.( 3.1l| ) produces the following recursive relation for the value 7„: 

^„ = ent/(") _ djn-l - " C'nln^s (3.12) 

completed by the initial condition 70 = 1. 

The localization transition in the n-particles model occurs at the values l3cP (j = 2, 3, 4, . . . , ?i) 
defined by the equation 

\ / 

where jj are the recursive solutions of Eq.( 3.12| ). 



3.3 Phase transition in the n— particle system for n ^ 1 

Now we are in position to give the answer to the question raised in Section "What should be 
the critical shape of the function /(n) (for n 3> 1) to have the simultaneous collective n-particle 
localization in the system?" For simplicity later on we suppose n to be even. The answer is very 
straightforward: The condition on the critical shape fcrin) providing the joint n-particle localization 
is defined by setting Uct'' = fain) in the equation ( 3.12| ) with -jj extracted from (l3.13| ). 



The approximate evaluation of the integral ( |3.13D for n ^ 1 gives us (with the exponential 
accuracy — see Appendix for details): 

7„ = exp ( — nlnn — — nln + O(lnn) 1 (3-14) 

We drop in Eq. (|3.13| ) all terms growing with n ^ 00 slower than exp(const n). 

Substituting ( p^ into ( |3H| ) it is easy to conclude that for n S> 1 the term 7„ has the maximal 
value in the sum (^.ll]) : 

where we kept the leading asymptotics (for n — > cxd) only. 

Estimating the sum ( pl|) with the same accuracy, we arrive at the following final conclusion: 

• If the function f{n) grows faster than fcvin) = |lnn then all n ^ 1 particles are localized 
simultaneously in the point well located at the origin x = of the halfline^^. 

Recall that [/^"^ is the depth of just the n-particle well. Thus the "integral" critical depth, 
Ucvin) = nllcr^ per the cluster of n particles (i.e. steps of different random walks at the point x = 0) 
has the following asymptotic behavior 

3 

Ucr{n) = -nlnn 

If the integral potential at the origin U{n) grows (with n 00) slower than Ucr{n) = |nlnn then 
only a finite fraction of random walks remains at the origin in the localized state. 
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4 Instead of discussion: Connection with random heteropolymer 
absorption problem 



The above mentioned problem has a direct apphcation to the old problem of the random copolymer 
absorption at the point well. Namely consider the random walk on Z"*" with the point well located 
at the origin and suppose that the well depth depends on the current "time" (i.e. the number of the 
random walk step) at the point x = 0. The partition function G^v of such process one can write as a 
sum over all random walk paths 



N 



{a;i,...,xjv_i}i=l 



j j — 1 I / ^ 



(4.1) 



where N is the normalization constant, g{\xj — Xj-i\) is the local transition probability and the 
potential Uj{xj) is the random variable of j. Suppose that Uj{xj) allows the following representation: 



Uj{xj) 



El + S2 for Xj = 
for Xj > 



where £1^2 are some positive constants and aj is the random variable defined as follows: 

+1 with the probability p 
— 1 with the probability 1 — p 



(4.2) 



Write (|4.l| ) in form of recursion relations: 

QN+lix) = @Nix + 1) + @Nix - 1) {X> 1) 

QN+i{x) = e^^QN{x + l) (x = l) 

e^(x) = (x<0) 

^ @N=o{x) = Sx,0 

Performing the shift x ^ x + 1 rewrite ( [4. 3D as follows 

G^+i(x) = A,@Nix) + (e^^ - 1) S^^iQNix + 1) (x > 1) 
6^(2;) = (x = 0) 

Qn=o{x) = Sx,l 



(4.3) 



(4.4) 



The partition function Q]\f{x) is the random variable which depends on the quenched random 
pattern of realisations of the potential Uj (j G [Ij-^]) at the origin. To find the reliable thermody- 
namic quantity we have to average the free energy In over the distribution of all random sequences 
{'7i,(T2,... jCTat}. The corresponding computations we realise in the frameworks of the replica ap- 
proach. 
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Averaging the ns power of the partition function {Q^{x)) = ^j^{xi, . . . we get: 

n 

[e^^ - 1J^6x„iA^^^„^^„^cc„^n{xi, ..,Xi-i,Xi + 1, Xj+i, .., + 
i=l 

' 2\ " 

(e^^ - \J2^^^^^^J^^Xl,..,i^,..,:i^i,..,x„^N{xl,■.,Xi-l,Xi + l,Xi+i, ..,Xj-i,Xj + l,Xj+i,..,Xn) + 



((e^^ - l)") . . . 4„,i$7v(xi + 1, + 1) 



(4.5) 



The solution of Eq.( |4.5| ) reads (compare to ( 3.11| )): 



<I>Ar(xi,X2, . . . = <I>iv(xi)Z^'^'^(x2, . . . ,a;„) + all permutations+ 

<I>iv(xi, X2)^^'^'^(x3, . . . , Xn) + all permutations+ 
(^N{xi,X2,x^)Z^^'^'^{xi, . . . , Xn) + all permutations+ 
••• + 

'^n{xi, ■ ■ ■ , Xn-i)Z].^'^'^{xn) + all permutations+ 
$Ar(xi, ...,Xn) 



(4.6) 



where the functions ^n{xi, . . . ,Xk) satisfy the master equations {k S [1,''^]) 

$Ar(xi, ..,Xj) = Aa;i,..a;^ ^>Ar(xi, ..,Xj) + . ,l7i^A'(a;i + 1, ■■,Xj + 1) {xi, ..Xj} > 1 

$Ar(at least one Xi>2 and all others Xj = 1) = 7^ i} G [li j] (4-7) 

and by 7^- we denote the averaged values of the Boltzmann weights: 

7. = (|(e^iv _ _ i)i + (1 _p) (e^2 _ (48) 

Using the results of the previous Sections one can conclude that the averaged moments of the 
quenched heteropolymer partition function, (Gat) , (0^) , • • • , (©jy) exhibit the singular behavior at 
the set of points being the solutions of following equations: 

p^'> = 7, + 1 = p (e^i - 1)^' + (1 - p) (e"2 _ ly + i (4.9) 



where = 2 (see Eq.{^)) and for j G [2, n] is given by Eq.(p^). Recah that we restrict 
ourselves with the case: ei > and £2 > 0. 

The "true critical point" of the localization transition in all averaged moments of the partition 
function 0^ (in the thermodynamic limit N 00) can be obtained using the following simple pro- 
cedure. We fix some arbitrary value £2 and find the minimal value e"(e2,j) among all solutions of 



Eq.(U) for j £ [l,n]. 



It is easy to check that for all 1 < j < n and any arbitrary choice of £2, the minimal value £1^ 
corresponds just to j = 1. It means that all moments of the random copolymer partition function, 0]y, 



13 



averaged over the quenched disorder in monomer types diverge at the same point as the "one-particle" 
part (i.e. "annealed" copolymer partition function), {Qn)- 

Of course, our consideration has an obvious crucial shortcoming connected with the fact that the 
repHca approach presented above does not ahow us to take properly the limit n ^ 0. Thus, the 
computations performed in Section ^ cannot be regarded as a proof of the conjecture that the phase 
transition points of copolymers with quenched and annealed chemical sequences coincide. However, 



the consistency of our investigation with other speculations on that subject [g, 17| at least gives hope 
that our conclusion is correct. 
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A Appendix 



Let us estimate the value of the integral /„ = ^ (see Eq.(p?7|)) for n ^> 1 

tt" Jo Jo 1 — cos ki . . . cos A;„ 



Changing the variables /cj = qn and expanding the nominator and denominator of the fraction of ( A.l ) 
up to the first non vanishing term, we get: 

Passing to the n-dimensional spherical coordinate system we arrive after simple algebra at the following 
expression: 



In the limit n ^ 1 we find the asymptotic expression of the function /„ with the exponential 
accuracy: 

(3 1 vre'^ \ 

— — nlnn + — nln + O(lnn) 1 (A. 4) 

This asymptotic expression has been used in the derivation of the equation ( 3.15| ). 
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Figure Captions 



Fig.l. The schematic representation of the multiparticle absorption model. 

Fig. 2. The choise of the lattice for the 2D-random walk in the first quater (xi > 0, X2 > 0) and the 
possible local moves of the walker. 

Fig. 3. The lattice (xi > 1, X2 > 1) with zero's boundary conditions obtained from the lattice shown 
in Fig.^ by the shift xi^2 a^i,2 + 1- 

Fig. 4. (a) The auxiliary 2D-lattice used in the computation of the nonmultiplicative contribution 



Wn{xi, X2) — see Eq.(2.17); (b) The same lattice as in (a) but with another boundary condition at the 
"main corner". 

Fig. 5. The phase diagram for the two-particle catalytic absorption model. 
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